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Abstract. A new approach is developed to understand stability of a population and further 
understanding of population momentum. These ideas can be generalized to populations in ecology 
and biology. 

1. Introduction 

One commonly prescribed approach for understanding the stability of system of dependent 
variables is that of Lyapunov. In a possible alternative approach - when variables in the system 
have momentum then that can trigger additional dynamics within the system causing the system 
to become unstable. In this study stability of population is dehned in terms of elements in the 
set of births and elements in the set of deaths. Even though the cardinality of the former set 
has become equal to the cardinality of the latter set, the momentum with which this equality has 
occurred determines the status of the population to remain at stable. Such arguments also works 
for the other population ecology problems. 

2. Population Stability Theory 

Suppose |PAr(to)| he the cardinality of the set of people, PAr(fo), representing population at global 
level at time to, where P/v(to) = {ui,U2, ■ ■ ■ ,un}, the elements ui,U2, ■ ■ ■ ,un represent individu¬ 
als in the population. Broadly speaking, the Lyapunov stability principles (see lyrni) suggests, 
|-P/v(^o)| is asymptotically stable at population size U, if ||PAr(T)| — P| < e (e > 0 ) at all T when¬ 
ever T > to- In some sense, |P;v(to)| attains the value U over the period of time. Lotka-Voltera’s 
predator and prey population models provide one of the classical and earliest stability analyses 
of population biology (see for example, |JDM] I and Lyapunov stability principles often assist in 
the analysis of such models. These models have equations that describe the dynamics of at least 
two interacting populations with parameters describing interactions and natural growth. Outside 
human population models and ecology models, stability also plays a very important role in under¬ 
standing epidemic spread |AR| . In this paper, we are interested in factors that cause dynamics in 
Pn and relate these factors with status of stability. A set of people Pm( 0 ) = {wmi, Mm2)''' > 
where Pm(^o) C Pjq{to), are responsible for increasing the population (reproduction) during the 
period {to, s) and contribute to Pn{s), the set of people at s (if they survive until the time s). The 
set Qmi('S — to) = represent removals (due to deaths) from PN(fo) during 

the time interval (fo,s). Let R0(s —fo) he the period reproductive rate (net) applied on Pm(^o) for 
the period (to, s), then the number of new population added during (to, s) is R,^(s — to) |-Pm(^o)| • 
Net reproduction rate at time to (or in a year to) is the average number of female children that 
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would be born to single women if she passes through age-specific fertility rates and age-specific 
mortality rates that are observed at to (for the year to). Since net reproductive rates are futuristic 
measures, we use period (annual) reproductive rates for computing period (annual) increase in 
population. Let C^iis — to) = {wi,W2,--- be the set of newly added population during 

{to, s) to the set P^ito). After allowing the dynamics during {to, s), the population at s will be 


-Pv(^o) U CAr^(s — to) — Qmi{s — to) 

( 2 . 1 ) 

Note that, Qmi(s— to) C PAr(to)UCAr^(s—to), because the set of elements ••• 

eliminated during the time period (to, s) are part of the set of elements 

{ui, U2, ■ ■ ■ ,un,wi , W2, ■ ■ ■ , wni} and the resulting elements surviving by the time s are repre¬ 
sented in equation fl 2 . 1 l) . The element ui in the set fl 2 . 1 l) may not be the same individual in the set 
Fjv(to)- Since we wanted to retain the notation that represents people living at each time point, 
so for ordering purpose, we have used the symbol Ui in the set fl2.1|] . 

Using Cantor-Bernstein-Schroeder theorem |PS| . |CjVj(s —fo)| = |Qmi(s —^o)| if |CAfi(s —fo)| 
< |QMi(s-fo)| and |C7Vi(s-fo)| > |QMi(s-fo)|- If |C7Vi(s-fo)| = |QMi(s-to)| then the 
natural growth of the population (in a closed situation) is zero and if this situation continues 
further over the time then the population could be termed as stationary. Assuming these two 
quantities are not same at to, the process of two quantities ICtvJ and |QmiI becoming equal could 
eventually happen due to several sub-processes. 

Case I; ICatjI > |QmiI at time to- We are interested in studying the conditions for the process 
|CniI —t IQmiI for some s > to. Two factors play a major role in determining the speed of this 
process, they are, compositions of the family of sets [{Pm(s)} : Vs > to] and [{R0(s)} : Vs > to]. 
Suppose |Pm(si)| > |Pm(’S 2)| > ••• > |PM(’Sr)| but the family of {|R0(s)|} does not follow 
any decreasing pattern for some to < < S2 < ... < st < s, then ICatjI ^ IQmiI by the 

time St- If R(^(si) > R(/,(s2) > ... > R0(sr) for to < si < S2 < ... < st < s such that 
|R0(sr — st-i) |PM('Sr-i)| — Qmi{st — st-i)| 0 for some sufficiently large T > to and suffi¬ 
ciently small |R(/,(st — sr-i)|, then ICatjI —)■ \Qmi\ by the time st- Note that in an ideal demo¬ 
graphic transition situation, both these quantities should decline over the period and the rate of 
decline of \Qmi\ is slower than the rate of decline in |CjVj| because ICatjI > \Qmi\ at time to. De¬ 
mographic transition theory, in simple terms, is all about, determinants, consequences and speed 
of declining of high rates of fertility and mortality to low levels of fertility and mortality rates. For 
introduction of this concept see Eni and for an update of recent works, see [in|. Above trend 
of |Pm(si)| > |Pm(s 2)| > ••• > |PM(sr)| (i-e. decline in people of reproductive ages over the 
time after to) happens when births continuously decrease for several years. Following the trend 
R0(si) > R(/,(s 2) > ... > R0(sr) will lead to decline in new born babies and this will indirectly 
result in decline in rate of growth of people who have reproductive potential. However the decline 
in |R<^(s)| for s > to is well explained by social and biological factors, which need not follow any 
pre-determined mathematical model. However the trend in |R<^(t)| for t < to can be explained 
using models by fitting parameters obtained from data. During the entire process the value of 
\Qmi \ after time to is assumed to be dynamic and decreases further. If a population continues to 
remain at this stage of replacement we call it a stable population. The cycle of births, population 
aging and deaths is a continuous process with discretely quantifiable factors. Due to improvement 
in medical sciences there could be some delay in deaths, but eventually the aged population has 
to be moved out of {Pat}, and consequently, population stability status can be broken with a 
continuous decline in {ICAril}. 

Case H: ICatJ = \Qmi\ at time to. It is important to ascertain whether this situation was 
immediately proceeded by case I or case H before determining the stability process. Suppose case 


u: ue PN{to) U Cn^s - to) 
and u ^ Qmi{s - to) 

= {Ul, M2, • • • , Wa+Ai-Mi} 
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b) Possibility of global stability even tiiough some sub-populations satisfy |C/v.|-|<?«.|=0 


Figure 2 . 1 . (a) The cycle of all the cases could follow one after another and the 
quantity at which equality of and Qmi occurs determines the duration of the 
case II. (b) Some of the sub-populations which are not satisfied the equality of 
and Qmi is compensated by the other sub-populations which are satisfying either 
CaTi > Qmi or CaTi < Qmi- 


II is immediately preceded by case I, then the rapidity and magnitude at which the difference 
between |CjVi| and |QmiI was shrunk prior to fo need to be quantified. Let us understand the 
contributing factors for the set Qm- At each f, there is a possibility that the elements from the 
sets Cat^, P/v — Catj — Pmi Pm are contributing to the set Qm- Due to high infant mortality rates, 
the contribution of into Qm is considered to be high, deaths of adults of reproductive ages. 
Pm, and all other individuals (including the aged), P/v — Cni — Pm, will be contributing to the set 
Qm- Case II could occur when |C/vJ and \Qmi \ are at higher values or at lower values. Equality at 
higher values possibly indicates, the number of deaths due to three factors mentioned here are high 
(including high old age deaths) and these are replaced by equal high number of births, i.e. |P<^| 
and \Pm\ are usually high to reproduce a high birth numbers. If equality at lower values of ICatjI 
and \Qmi I occurs after phase of case I then the chance of P/v remaining in stable position is higher. 
Suppose elements of P/v are arbitrarily divided into fc—independent and non-empty subsets, A(I), 
A{ 2 ), ■ - - ,A{k) such that IPa^H/^ |y 4 (s)|ds. Let P be the family of all the sets A(s) such that 

U(A(s)) = P/v. Members of P are disjoint. Suppose ^ ^ be an arbitrary size of k* of subset 


of P are satisfying the case II and P 


P 

k* 


are not satisfying at time to and t > to, then we 


are not sure of total population also attains stability by Theorem [T] 
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Theorem 1. Suppose each of the member of 
and F — 


F 

k* 


F \ 

j are not satisfying the condition \ \ 

always leads Pn to stability. 


is satisfying the condition | CVJ = \Qmi\ 
= \Qmi\ o,t time t > t^, then this does not 


Proof. Note that F has collection of fc—sets. Suppose a collection C divides into fc—components 

of subpopulations {CAr^(l), C7Vi(2), ■ • • ,CNj^{k)} such that |C7vJ=|CAri(s)| ds, where C]s[^{s) is 
the subset in C and a collection Q divides Qmi into k— components of subpopulations 

{Qmi(1),Qmi(2), • • • ,QMi{k)} such that 

IQmiI = /I |C'mi(s)| ds, where Qmi{s) is the subset in Q. 

By hypothesis, |CArj(s*)| = |Qmi(s*)| for s* G {!*, 2*, • • • , k*} at each time t > to until, say, tr- 
The order between k* and k — k* could be one of the following: 2 k* < k, 2 k* > k, k* = ^. Suppose 
Cni F C and Qmi C Q with 


c;. = {c;.(i).c;,(2),...,c;.w} 

Qmi = Qmi(2)) • ■ ■ ) 

for same above arbitrary combination of /c*—components and rest of the k — k* components are 
satisfying |C^* (s**)| — |Qmi(s**)| 7^ 0 for all s** = 1,2, ■■■, k — k*. We obtain unstable integral 
over all k — k* components to ascertain the magnitude of unstability. 



The stable integral for this situation is 

(2.3) f [|CJ,,(s-)|-|Q*„,(s*)|]*- 

To check the unstable and stable points over the time period {to, tx), one can compute following 
integrals; 

(2.4) 



(2,5) F f [|CJ,.(s*)| - [Oy (s-)|] ds’du 

Jto Jl 

For each of the k — k* component, the values of (s**)| — can be either positive 

or negative. If at time F, for all s** = 1,2, ■ ■ ■ , k — k*, the values of 

are positive (or negative) then the eq. fl2.2p will take a positive (or negative) quantity and the 
population at time to is not stable. If such a situation continues for all tx > t^, then the integral 
in eq. fl2.4p would never become zero and the population remains unstable in the entire period 
(fo,^r)- However, for some of the s**, if the quantity (s**)| — is positive and for 

other s**, if the quantity |C')^*(s**)| — is negative such that eq. fl2.2p is zero at each 

of the time points for the period (to)^r) then the population remains stable during this period 
(because by hypothesis the eq. fl2.5p is zero). □ 
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Case III. ICatJ < |QmiI at time to- Global occurrence of this case at lower values of ICatJ and 
\Qmi\ indicates that the is declining and also is in unstable mode. R^p has been very low 
consistently for the period t < to and the supply to the set Pm has diminished over a period in the 
past. All the subsets of and Qmi might not be stable in case III, but by similar arguments 
of the Theorem [H global population behavior nullihes some of the local population and case III is 
still satished globally. 

All three cases would be repeated one following another. Most countries are currently facing 
case I with varying distance between ICatJ and |QmJ • 

3. Replacement Metric 

We introduce a metric, which we call a replacement metric, with a space, Mr as follows: 

Definition 2. (Replacement Metric). 

Let Ai = min {||CAri(s)| - |Qmi(s)|| : s > R} and 

A 2 = max{||CATi(s)| - |Qmi(s)|| : s > to}- Let Mr = [^ 1 ,^ 2 ] C M+ and M = {||Cati(s)| - |Qmi(s)|| 
: s > to} with the metric dM{x,y) = . We can verify that {M^du) is a metric space with 

dM '■ {M X M) —)■ Mr and non-empty set M. 

The metric M, in the dehnition 1 is bounded, because dM {x,y) < k for k > 0. 

Definition 3. Suppose ||C7Vi(si)| - |Qmi(si)|| = /i, ||CAri(s 2 )| -\Qm 2 {s 2 )\\ = /2 and so on 
for Si < S 2 < ... . Then we say population is stable if fsj, —>■ 0 for sufficiently large T and 
^ |Civi(st)| = ^ |Cvi(st)| = 0. 

4. Conclusions 

We can prove that the value at which the population remains stable is variable, i.e. the value at 
which the population becomes unstable by deviating from case II could be different from the value 
(at a future point in time) population becomes stable when it converges to case II. Replacement 
metrics (see dehnition [2]) are helpful in seeing this argument and such analysis is not possible by 
Lotka-Voltera or Lyupunov methods. Due to population momentum, there will be an increase in 
the population even though the reproduction rate of the population becomes below the replacement 
level. Population stability will always attain a local stable points before diverging and again 
converging at a local stable point. The duration of a local stable point depends on the density of 
the population and resources available for the population. 
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